CONVERGENCE OF THE EIGENVALUE DENSITY FOR 
/3-LAGUERRE ENSEMBLES ON SHORT SCALES 
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Abstract. In this note, we prove that the normalized trace of the resolvent of the /3- 
Laguerre ensemble eigenvalues is close to the Stieltjes transform of the Marchenko-Pastur 
(MP) distribution with very high probability, for values of the imaginary part greater 
than m"^"''^. As an immediate corollary, we obtain convergence of the one-point density 
to the MP law on short scales. The proof serves to illustrate some simplifications of the 
method introduced in our previous work 20^ to prove a local semi-circle law for Gaussian 
/3-ensembles. 



1. Introduction 

Consider anmxn matrix X, whose entries are i.i.d. complex Gaussian random variables 
with mean and variance EjXjjp = 1. The m x m matrix H = XX*, the star * denoting 
the conjugate transpose, is a Wishart matrix [23]. It is a classical result of random matrix 
theory that the distribution of the eigenvalues Ai, . . . , of H is given by the density /2(A) 
on M™: 

m 

MX) = n - ^/ • n Ar-e-Z^Sr^i 

i<j i=l 

Here Z2^m is a normalization factor. Suppose the limit m/n — )• d exists for < d < 1. 
Then the empirical distribution 

F2,m{x) = — • tt{l < « < : Ai < x} 
m 

of the eigenvalues of the rescaled matrix H/m converges to the Marchenko-Pastur distri- 
bution, with density 

pMP,d{x) = ^^x/(A+ - x)(x - A_)l[A_,A+](a;), 

where 

A± = (l±\/d)2. 

In fact, V.A. Marchenko and L.A. Pastur [T2] showed that weak convergence of the eigen- 
value distribution also holds when the entries of X are i.i.d. but not necessarily Gaussian. 

A. Edelman and I. Dumitriu [7J have introduced an infinite family of tridiagonal random 
matrix models, termed /3-Laguerre matrices, which generalize the Wishart model, and have 
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the explicit eigenvalue density 

m „ 

(1) W(A) = Z.i^HlX. - A/ . n Ar-^-^-^^-^e-^^r.. 

i<j 1=1 

for < /3 < oo. Here a > (3{m — l)/2 is a real parameter. Zp^rn is another normaliza- 
tion factor. The (Gaussian) Wishart eigenvalue distribution corresponds to the Laguerre 
ensemble with (3 = 2 and a = /3n/2. Almost-sure convergence of the eigenvalue density to 
the Marchenko-Pastur distribution with parameter 

d = lim — — < 1 
2a 

was established for the eigenvalue distributions of the /3-Laguerre matrix models by the 
moment method in [6]; that is, for any a < b: 

Ff5,n{b) - Fi3^ri{a) / PA/P,d(s) ds, 



almost surely as n — )• oo. Fp^n{^) is the eigenvalue distribution function, defined as in the 
case (3 = 2 above. 

The object of the present note is to extend these results on convergence of the eigenvalue 
distribution of general /3-Laguerre ensembles to short intervals [om, bm] such that bm — am = 
0{m~'^^''), for e > arbitrary. 

The main result is the following: 

Theorem 1. Let /3 > 0, and < d < 1. Set a = m(3/2d. Let 5,K.,e > be positive 
parameters, and E G (A_ + k,A_|_ — k). Denote by Nj the number of eigenvalues of the 
(3-Laguerre ensemble of size m with parameter a in the interval L. For any k > 0, there 
exists a constant Cs^K,e,k such that: 

P — - / PMP,d{x) '^^ 



The proof of Theorem [T] is obtained by combining Theorem [2] with Corollary [3] in Section 
m Following recent work of L. Erdos, B. Schlein, H.T. Yau and collaborators in the case 
of Wigner matrices (see for example [9], [10], [11], [12], [13], [H]), our method is based on 
the study of the resolvent matrix of the Edelman-Dumitriu tridiagonal models. 

Let M/3^rf denote a normalized tridiagonal /3-Laguerre matrix (see Section [2] for details 
on notation). The imaginary part of the trace of the resolvent 

spA^) = ^ tr(M^,d - z)~^. 

provides an approximation for the eigenvalue distribution on scales comparable to the 
distance between z and the spectrum of Mp^^- In [22]) we showed how a resolvent expansion, 
together with an iterative argument based on the Schur complement identity could be used 
to derive a local version of the semi-circle law for the Gaussian /3-ensembles. Tridiagonal 
models for these eigenvalue distributions appeared in [7] . In Proposition 12. 1^ we use a 
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resolvent expansion to show that sp^d{z) is close to the Stieltjes transform of PMP,d for 
values of z away from the spectral edges and with 9z > m^^^'^^''. This implies that 
Theorem [1] holds for intervals / of size |/| > m~^/^+'^. The argument in the present 
work is substantially simpler than the corresponding one in |20j, where we proved a local 
convergence result on the scale m~^^'^~^'' using a resolvent expansion and asymptotics for 
Hermite polynomials derived by the Riemann-Hilbert method. Instead of attempting to 
exploit the cancellation due to oscillation of the normalized Laguerre polynomials, we use 
a general off-diagonal resolvent estimate (see Lemma [22]) • The computation of the limit of 
the normalized resolvent trace for the deterministic matrix corresponding to /3 = oo using 
Riemann-Hilbert asymptotics in [20] has been replaced by a less involved derivation, see 
Section [31 

The iterative argument leading to Theorem [2] in Section [4] is similar to the one in [20] . 
Note that Theorem[2]is deduced from Propositions 14.31 and 14. 41 without reference to Propo- 
sition [2TTJ In contrast, the local result on the intermediate scale m"^/^"*""^ in [20] was used 
as an input for an inductive argument to reach the scale m~^"'"'Q. Although it can be 
entirely replaced by the iteration in Section [4] we have chosen to present the argument 
for Proposition 12.11 because it provides an elementary alternative to the moment method 
or the Schur complement approach to proving the Marchenko-Pastur law. The proof of 
Proposition 12.11 does not depend on the specific properties of the /3-Laguerre ensembles 
other than concentration of the entries around their mean. It can be applied to fairly 
general tridiagonal models with independent entries to prove convergence of the eigenvalue 
distribution, and convergence up to some intermediate scale depending on the magnitude 
of the entries. We give examples of such extensions in Section [5] 

We end this section with some references to previous literature. Local versions of the 
Marchenko-Pastur law for the eigenvalue distribution of XX* when the entries of X are 
independent but not necessarily Gaussian have appeared in [1], [12], [21]. The first paper 
deals with the hard edge of the spectrum in case d = 1. The use of a perturbative expansion 
around deterministic matrices associated to Laguerre polynomials already appears in [8], 
where the authors study fluctuations of the spectrum in the large /3 limit. In [19j, I. Popescu 
proves convergence and Gaussian fluctuations for the moments of tridiagonal matrices 
under a scaling assumption for the moments of the entries. In contrast to the result in [20] 
at the time of publication. Theorem [1] appears to be new for general f3. 



2. Tridiagonal models and resolvent expansion 

We recall the main result from the work of A. Edelman and I. Dumitriu [7], which will 
be our starting point. Let m G N, /3 > and choose < d < 1. Let a = mj3/ (2d). Consider 



Modulo minor changes, the inductive argument in [20] can also be used to obtain the semi-circle law 
down to scale 72"^+". 
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the bi-diagonal matrix 



Br- 



( X2a 

Xp{m-1) X2a-(5 

V 



X/3 X2a-/3(m-l) / 



In the equation above, the symbol Xri f > ^ represents a random variable with chi distri- 
bution with r degrees of freedom, defined by the probability density function 



r-l - — 

X e 2 . 



r(r/2) 

r denotes the Gamma function, defined for x > by 



r(x) 



t 



x-l. 



dt. 



The random variables appearing in the matrix i?^ are all independent. By [7], the eigen- 
value distribution of the real tridiagonal matrix = BpB^p is given by the density 
function /^ ^ in ([T|. We will be concerned with the scaled /3-Laguerre matrix model, de- 
fined by 

We will determine the behavior of the normalized trace of the resolvent, S/3,d(-z) in a neigh- 
borhood of the real axis in the upper half-plane C''" = {^z > 0}: 



1 



Si3,d{z) = — tr {Mis^d - z) 



-1 



m 



The imaginary part 3's^^d(z) is the Poisson integral of the eigenvalue distribution. As such, 
it is an approximation to the empirical eigenvalue distribution of Mp^d at scale To prove 
Theorem [U it will suffice to show that, with probability no less than 1 — Cs^K,e,kfn~^ , we 
have 



where 



SMP,d{z) 



z + 1 



X — z 



2dz 



d V{z-X.){z-X+) 
2dz 



is the Stieltjes transform of the Marchenko-Pastur distribution with parameter d, (see [13 
Lemma B.l or [T3] 7.1, as well as Corollary [3] below) . 
We begin by writing 



(2) 



I3,d 



ooA 



+ A, 
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where M^o^d is the matrix 



d 
m 



( mid ^mld\/ra — 1 

^mld\/m. — 1 m/d + m — 2 ^Jm,|d — l\/m — 2 



\/m/d — m + 2\/2 m/ d — m + A ^ mj d — m ^ 1\/T 

\ \Jmjd — m + 1\/T m/d — m + 2 ) 

The entries of the tridiagonal matrix A are 

/3m /3m 



. ,X2a-(i-l)X/3(m-j) " j - - j 

■^■^ m 
For large m, all the entries of A are simultaneously small in magnitude, with overwhelming 
probability: 

(3) \Ajk\ < m-V2+^ l<j,k<m. 

for any c > 0. Here and below, we will say an event E = E{rn) holds with overwhelming 
probability if, for each k, there is a constant Ck such that 

This follows readily from the definitions, and properties of the Xr cind Xr distributions, 
which are concentrated around their mean 

with exponential tails. 

2.1. The matrix M^o^d £ind generalized Laguerre polynomials. The spectral theory 

of the symmetric matrix Moo,d can be described explicitly in terms of Laguerre polyno- 
mials. For a > —1, the generalized Laguerre polynomials L^, A: = 0, 1, . . . are orthogonal 
polynomials with respect to the measure 

Wa{x) dx = x°'e~^ dx 

on the positive real axis M"^ = [0, oo). is normalized such as to have L^(M+, uiadx)-norm 
1: 

Lk{x)L'^{x)wa{x) dx = Sj^k- 



I 



The eigenvalues of Moo,d are the normalized zeros Aj = d ■ li/m of the mth generalized 
Laguerre polynomial with parameter 

,1 

a = m ■ - — 1 
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There is a complete set of correspond eigenvectors vi of the form 



(4) 



L 



m(l/d-l)+l 
m-2 



\ ^0 



for < i < m — 1 . We denote by Ui the normahzed eigenvectors 



Ui 



Vi 



Vi 



This follows from the recurrence relations for the classical Laguerre polynomials L^, nor- 
malized to have leading coefficient 

The polynomials satisfy (see [21 22.7.30-22.7.31]): 

(5) 
(6) 



-[{k + a + l)Ll{x) 



{k + l)Lt^^{x)]. 



2.2. The resolvent expansion. In this section we show that the normalized trace of the 
resolvent is well approximated by the Stieltjes transform of the Marchenko-Pastur density 
for 9z > m""^/^^'^, e > 0. Specifically, we have the following 

Proposition 2.1. Let e,/t > 0. For any z such that 9z > m"^/^"*"^ and 

A- + K < 5Rz < A4_ — K, 



we have 



\sp{z)- s^{z)\ =0,,^{m-'l^) 



with overwhelming probability. 

The choice of the scale m~^^'^~^^ in Proposition 12.11 is largely arbritrary, given that the 
inductive argument of Section fallows one to reach the optimal scale m"^"*"*^. The exponent 
I — e was chosen because it represents a scale accessible without detailed information on 
the size of the entries of the resolvent. Using the estimates for Laguerre polynomials in 
|16j . and the method of |15] to obtain a lower bound for \\vi\\, one can prove sharper bounds 
than ^ for |(Moo,d - z)~j^\ when \i - j\ < m^/^. The expansion ([8|) can then be summed 

for values of Qz smaller than m~^/^, but still much larger than m"^"*"^. We will not pursue 
this here. 

The proof of Proposition 12.11 proceeds by a resolvent expansion comparing the trace 
of the resolvent of M^^^ to that of the deterministic matrix Moo,d. The precise spectral 
information available for the matrix Moo,d allows us to calculate the large m limit of the 
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normalized trace Soo{z), and to control the resolvent expansion ([8]). Starting from ([2]), 
write: 

(M^,d - z)-^ = (Moo,d - z)-^ - (Moo,d - z)-^A{Mfs,d - z)-\ 
Upon iteration, this yields 

(7) {Mp^d - z)-^ = {M^,d - z)-' 

i-i 

+ {-{Mo.,d - zy'Af {M^^d - z)-^ + (-(Moo,d - z)-^A)\Mp^d - z)-\ 

k=l 

Taking traces and normalizing: 

(8) sp^d{z) = s^{z) + - 5^ tr (-(M„o,d - z)-^Af {M^^a - z)-^ 

k=l 

+ 1 tr(-(Moo,d - z)-iA)'(M^,d - z)-\ 
m 

By (fT7|) . we can replace the first term on the right with SMP,d{z), introducing an Oi^{m~^) 
error. It suffices to show that the remaining terms also decay at least at an 0{m~'') rate 
as m — 7- oo. 

For i < i, j < m, we let 

Rij=Rij{z) := {M^,d-z)-j\ 
Expand the A;th term in of the sum in ([8]) to find: 

l<il,...ife+l<rra 

where = i;, + 1 or — 1. Taking absolute values and using ([3]), this is bounded by 

(9) ^-(/c+2)(l/2-c) ^ |i?,^,,||/?,,,,3|...|i?,,^^,J, 

l<ii,...ifc+i<m 

with overwhelming probability. 

The next lemma gives an estimate for the entries Rij of the resolvent which is efficient 
when i and j are widely separated. 

Lemma 2.2. For any 1 < i, j < m, we have 

(10) \R,j{z)\ < C{Qz)-' ■ exp(-Q(9z) • \i - 

Several results on exponential decay of the resolvent entries could be used to obtain 
Lemma [221 We will use a slight variant of the estimate of Combes-Thomas type developed 
by M. Aizenman in the context of localization for discrete random Schrodinger operators: 
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Lemma 2.3 (M. Aizenman, p], Lemma ILl). Let T be a countable set with a metric 
d : r X r — )• [0,00), and let H be a self-adjoint operator on i'^{T) whose off-diagonal 
elements are exponentially summable: 



fill 



5„ = sup V|i^(x,y)|e"'^(^'^) < 00. 

X 



Then, for energies not in the spectrum of H , with 

A = dist(z, o-(i?)) > 0, 

we have 



aA 



d{x,y) 



Note that the hypothesis (jlip differs from the one in Lemma ILl in [T], but the same 
proof works also for Lemma 12.31 as stated. 

To bound (j9]), we can assume that ij = i'j. Indeed, the bound (llOp only changes by a 
constant factor when the index i is changed by a unit. Consider the sum over in ([9|): 



n=i 



We split this sum into three regions: 



+ 5^ + 5^ ) l-^n«2 1 l-^ifc+iil I 

|ii-j2|>ml/4 |ii-i2|<ml/4_|j^_j^^j|>ml/4 |j^_j2|<ml/4^|jj_j^^-^|<ml/4 , 



By pOj) . the first two sums are bounded by 

Cm^^"^ exp(— m*^). 
We are reduced to considering the sum 

(12) l-^nj2 1 I- 

|n-j2|<rrti/4 

The second factor in each summand of (1121) is estimated as follows: 



(13) l-^iiifc+il - X] 



Ui{ii)ui{ik+i) 



i=l 



Xi- z 



< 



E 



9\ 1/2 / m 9\ 
' \Ui{ik+l)\ \ 



Ai - z\ 



|Aj - z\ 



The first factor on the right side of the inequality ()13p is no greater than m^^^ "^Z^, since 
lluill = 1 and 
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In summary, (|T2|) is bounded by 

l/8-e/2 



m 



We proceed to sum over ^2 in ([9]): 

We can now repeatedly sum over is, . . . ,ik-i, using (jlOp at every step. Each summation 
results in an additional factor of m^/^"*^. Thus ([9]) is bounded by 

(14) c'=-W/«-^/2.n^(^-i)/-(^-i)^5:K;M^ . 



Using the Cauchy-Schwarz inequality as in (fT3]) . we have 
By ([HD and CoroharyEl 

yy\!^p±i)l = y-l—<cmiogm. 

\\i-z\ ^\\i-z\~ 

Inserting the previous two inequalities into (fT^ . we find that the sum ([9]) is no greater 
than 

C^logm•m(^+2)/2-'=^ 

Choosing c smaller than e/2 in ([3|), this last quantity is 0{m~^'^l'^) with overwhelming 
probability. As for the final term in ([8]), it is bounded by 

^-(.+2)(i/2-c) ^ ,,3| • • • ^^,J • 

l<zi,...ifc<'m 

smcG 

|i?g(z)| = |(M;,,,-%|<mV4- 

for each i and j, provided > m^/^""^. Performing each of the k sums using (jlOp as 
previously, we find that the sum is bounded by 

with overwhelming probability. Letting /c in ([8]) be larger than 8/e, we obtain 
(15) Sfi{z) = s^{z)^0(rn-'l''). 

with overwhelming probability. Combined with the approximation (jl7|) . the relation (jl5p 
implies Proposition 12. 1[ 
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3. Convergence for /3 = oo 
In this section, we identify the Hmit of the resolvent 

'-'oo \ ^ ) — / \ — ~ 

^fr{^j-z m L?^{mz/d) 

for z as close as m~^^'' to the limiting spectrum, but away from its edges. In [20], the 
corresponding quantity for Hermite polynomials was estimated using asymptotics for these 
orthogonal polynomials in the complex plane. Here, we use a differential equation for Soo(-z) 
derived from the ODE satisfied by Laguerre polynomials. 

+ (« + !- x)(L° )'(x) + mL^ix) = 0, 

with a = m{l/d — 1). After rescaling, we obtain an equation for Sis{z) (see [8^ Section 
2.2.2]), 

1 dz 

(16) zd{soo{z)f + Soo{z) (d - 1 + z) + 1 + —Soo{z) H Sooiz)' = 0. 

m m 

To solve the equation approximately, we treat the final two terms 

e = —Sooiz) + —Soo{z)' 
m m 

as error terms, and use the rough estimates 

\soo{z)\ < m~\ \s'oo{z)\ < m^~\ 
so that £ = 0{m~^). For z in the region 

{Qz > 0,\z - A±| > K, \z\ < 4}, 

we find the solutions 



I , , z + d-l±'x/(z- X+)(z - A_) ^ , 



The solution s satisfies 9s > for 9z > 0, and we conclude that the right side of the 
equation is an approximation for S/j^^: 

(17) Soo{z) = smpAz) + 0«;(m"'). 

We expect that the error term can be replaced by 0(l/m), but ()17p is sufficient for our 
purposes. 

4. Inductive argument 
In this section, we improve the convergence of the density on short scales from the level 



m 



-i/4+e ^YiQ optimal level of m ^"'"'^ by an inductive argument. 
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Theorem 2. Let s{z) = sp^d{^) he the Stieltjes transform of the measure induced by the 
eigenvalues of the normalized j3-Laguerre ensemble matrix M^^. Let SMP.diz) he the Stielt- 
jes transform of the Marchenko-Pastur distribution. Then, with overwhelming probability, 

(18) sup \s{z) - SMP,di^)\ = o(l) 

where the domain D is defined as 

Z)e,K := {z -.Qz > m"^+^ A- + K < JRz < A+ - K, \z\ < 4}. 

To prove Theorem [21 we need three facts about the tridiagonal models and Stieltjes 
transforms. The first can be found in [6j: 

Proposition 4.1. One can diagonalize 

Mp,d = QAQ* 

such that the first row of Q is independent of A and consists of independent entries with 
Xp distribution, normalized to unit norm. 

The next corollary establishes the link between control of Stieltjes transform and control 
of the eigenvalue density. See for example [13j : 

Corollary 3. Let z = E + ir] and m~^^^ < < \E, < E <1G and e, r > 0. Suppose 
that one has the bound 

\sp,d,m{z) - SMP,d{z)\ < T 

with overwhelming probability for all such z. Then for any interval I in [A_ + e, A+ — e] 
with \L\ > r], one has 

Ni - n j pMP,d{y) dy < CeTn\L\ 
with overwhelming probability, where Nj denotes the number of eigenvalues in I. 
Finally, we shall need the following standard result. A proof can be found in |20j : 



Lemma 4.2. Suppose the Marchenko-Pastur distribution holds at level m" for some —1 < 
a < 0, that is, for any c > 0, 

\s{z) - SMP,d{z)\ < C 

for sufficiently large m, with overwhelming probability. 
Then we have 

(19) 1 ^ C{^z)-^m- log m 

for any z such that m^^ < 9z < m°. 

The proof of Theorem [2] is a combination of the following propositions: 
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Proposition 4.3. Suppose the Marchenko-Pastur law holds at level m"' for some —1 < 
a < 0, and that we have 

(20) \R^^f{z)-SMpA^)\=0{l) 

for > m° with overwhelming probability. Then, for any 6 > 0, 

(21) \R^,f{z) - smpA^)\ = o{l) 

holds for z such that 9z > m^""^)/^^'^ with overwhelming probability. 
For a given sequence X„ of random variables, we write 

Xn = o(l) 

with overwhelming probability if, for every c > 0, the event {|^n| < c} has overwhelming 
probability. 

Proposition 4.4. Suppose \Rif'{z) — SMP,d{z)\ = o(l) for z such that > m"" for some 
— 1 < a < 0, then we have an improved Marchenko-Pastur law, that is: 



(22) 



— yZ "r~ SMP,d{z) 

m ^-^ Xj — z 



0(1) 



for z such that > m^"" i)/2+<5 for any 6 > with overwhelming probability. 

Proof of TheoremU^ The condition (j20p holds trivially if Qz > 1. The result follows by 
repeatedly applying Propositions 14.31 and 14.41 □ 

What remains is the proof of the two propositions above: 

Proof of Proposition By Schur's complement, we have the following relation 



(23) R^fiz) = ^ . 

2 22 -Rii 

«ii - ^ - aiia2i i_^_„2^ jj^^ 
Here an denotes the normalized (1, l)-entry of the bidiagonal matrix i?^, 02,1 the nor- 
malized (2, l)-entry of Bp, R is the resolvent of := -^BpBi^ where Bp is formed by 

2 1(7' <^ 1 cri- V> 1 T-n 1 -i- t-1 ci 

/3m ' 



removing the first row and column of Bp. We remark that a\i is distributed as ^X2a and 

all distributed like ^X/3(m~i)- Lastly, an, 021 and R are independent. 

By the argument in Section [Sj we need only show that R^f satisfies the approximate 
functional equality 



(24) <(^) + . ^ ^ ' , = o{l), 

d — 1 + z + zd ■ R^{ [z) 



with overwhelming probability for z such that G [A_ + k, A+ — k] and > m^"^^)/^^'^, 
with 6 > arbitrary. We first note that, due to the restriction on ?R.z (in particular, 5?z > k 
since A_ > 0), equation (p^ is equivalent to: 

(25) zd ■ {Rii{z)f + {d-l + z)-Rii{z) + l = 0(1) 
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where from now on we suppress the superscript /3, d. 
Rewriting equation (j23|) . we have 



(26) 2021 • ^11^11 + zRu + aii-Rii - al^Rii + 1 = 0. 
We can further rewrite the above as 

Tn — 1 ' — - — ]^ ^ — - 

(27) zd RuRu + zRii + d Ru - Ru + l 

m m 

+ [z ■ rTiRu + RTi) (ail - ^^^^) - («n " 1) " = 0- 
It suffices, therefore, to show that 

7Tl — 1 - — - 

(28) Ru = Rn + 0(1) 

m 

and 

(29) (z • rTiRii + rTi) (^21 - ^^^^) - («n " 1) " = '^(l) 

for the set of z we are interested in. 

We first note that a^i has distribution ^X^(m_i) ^iid so 



a2i - d- 



m 



with overwhelming probability for all 6 > 0. Similarly, Cii has the distribution of -^X^rn/d 
and so 

|a?i - 1| < C7m-i/2+5 

with overwhelming probability for all 5 > 0. The assumption of the proposition implies 
that for all w with > m'^, we have 

\Rii(w)\ <C, 

for some constant C. By ([28]) . which will be proved independently below, we also have 

\Rri{w)\ <c. 

The inequality 

|i?ii(^)| < '^\Rii{w)l 



for S^z < Qtu, now implies that approximation (I29p holds with overwhelming probability 
whenever 9z > 77i("~i)/2+'5. 

We turn our attention to (128p . Firstly, we write 
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and similarly 



where (gi, . . . qm-\) is the first row of the eigenvectors for M^, Aj are the eigenvalues and 
(gi, . . . , qra) is the first row of eigenvectors for Mg, and Aj the eigenvalues. 
We write 



Rii • Rii = > Eo 



m — 1 
m 



^ Xj — z ^ ^ Xj 

3 3 



Eql m — 1 _ X - g? 

3 ■' 3 1 

where Ex denotes the expectation with respect to the random variables X. 

By Proposition 14. H q and Xj are independent and so are q and Xj. We write the first 
term in (I30p as 

.2 .2 2_1_ 

^1 m 



^ A,- - z ^ ^ A,- - z ^ Xi-z 

3 3 3 

Since q and Aj are independent, we can condition on Xj and apply Proposition 14.51 be- 
low to the sum to conclude that the right side of the previous equation is bounded with 
overwhelming probability by 

1/2 

By Lemma 14.21 the latter is bounded by C{'^z)~^m~^^'^^"'^'^~^'^logm for some constant C 
and any c > 0. The fifth and sixth terms of equation pop are similarly bounded. To deal 
with the middle two terms, we use the interlacing property of the eigenvalues of a matrix 
and its minor, and then the interlacing property of the eigenvalues of a matrix and a rank 
1 perturbation, to obtain that these terms are bounded by 

1 

m^z 

We end the proof of Proposition 14.31 with the following simple variant of McDiarmid's 
inequality, a proof of which can be found in the appendix of |20j : 

Proposition 4.5. Let Xi, . . . ,X„ be independent subgaussian random variables. Let Q C 
M" be such that 
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with overwhelming probabibility. Let F be a real function of n real variables such that if 
X'Y^ • • • 5 oCrp^^ Xj^ ^ J then 

I , . . . , Xn) . . . , Xi—i, Xi, ^i+l ) • • • 1 Xn)\ ^ Cj 

for all 1 < i < n, and outside of 17, F is bounded by a polynomial in n. Then for any 
A > 0, one has 

P(|F(X) -E(F(X))| > \a) < Cexp(-cA2) 
for some absolute constants C,c > 0, and a = c?. 

□ 

Proof of Proposition \4-4\ The proof is similar to that of Proposition 14.31 By the assump- 
tion, it suffices to establish that 



m ^—^ Xj — z 



0(1) 



with overwhelming probability. The difference inside the absolute value sign is ■ 
The statement now follows by another concentration argument and Lemma |4.2[ □ 



5. Some extensions 

The resolvent expansion in Section [2] allows for a quick proof of convergence of the 
eigenvalue distribution for tridiagonal matrices with independent entries "close to" a Jacobi 
matrix whose limiting spectral density is known. Consider a sequence of random tridiagonal 
matrices given by 

/ On bn-i \ 

bn-l On-l bn-2 



(31) 



\ 



ai bo 
bo ao J 



where a,-, bj, 1 < j < n. Suppose that 



(32) 
(33) 



0(n~ 
= 0(n 



uniformly in 1 < j < n with overwhelming probabilty, for some a > 0, where 
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are such that j~"bj, j~°^aj are of order 1. If the deterministic tridiagonal matrix has 
a hmiting spectral density fi, that is, if the spectral measures of n""^ • converge 
weakly to some measure /i, then 



- tr (n'^EAn 
n 



1 



z-\ 



^(dA) + o(l) 



for z in any compact subset D of {9z > 0}. Note that by (j32|) . (j33|) . the quantity Sa 
in Lemma 12.31 is bounded by n''/'^ with overwhelming probability ii H = n~"An. Using 
Lemma 12.21 as in the proof of Proposition 12.11 to perform a one-step resolvent expansion, 
we have, for z £ D: 



(34) 



1 



n 



tr (n-°A„ 



1 



z-X 



/i(dA = o(l) 



with overwhelming probability. Thus the eigenvalue distribution of n^°'An almost surely 
converges weakly to fj,. 

As an application, consider the following example from [19]. Consider the Jacobi matrix 

/ (n-l)° \ 
(n-1)" (n-2)" 

2)" 

2" 

2" 1 
1 



in 



V 



The moments of the rescaled matrix n "Jn converge to those of the Nevai-Ullman distri- 
bution [22], [T71 Section 5]: 

^■tiin-^Jn)" 



n 



fj-k 

Ua{dx) 

By a standard density argument 



-2,2] (2;) 

' air 



^ r-2 ^-l+l/a 



At. 



x\l2 \/4^ 

-tr(n-"J„-z)-i = / -i— z.„(dA) + o(l). 

CLj — {"^n^jj ~l~ 

where Xj, 1 < j < n — 1 and Yj, 1 < j < n, are mean zero random variables whose kth 
moment is bounded uniformly in j for all k, then for any /3 < a, (j34p holds for the matrix 



If we let 
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(35) 



n 



E 



1 







sup n'^E|an| < oo. 



for some e > and k = 1,2, ... , \^~\ . Then for any c > 0, we have 

2 



> cn 



Oc(l), 



uniformly in j, and a similar bound for aj. The convergence (j34p with fj, = holds in 
probability. Popescu assumes ()35p holds for all k with e = and obtains convergence 
of the moments of n~°'An as well as almost sure convergence to the limiting distribution. 
Note that [19] contains more general results that apply also to cases not easily accessible 
by our method. 

Given more precise information on the size of the entries of An or the rate of convergence 
of the Stieltjes transform for the deterministic model, one can improve on the above. To 
give a simple example, we introduce a "positive temperature" version of the Jacobi matrices 
associated to the orthogonal polynomials Pn{x), n > 0, for the measures 



(36) 

Define the density 



-''"dx. 



gn,m{x) — l[0,oo) T 



Let the random variable Xj have distribution given by gj,m, with all the Xj independent. 
We have 



EXi 



r i+1 



2-l/2m / j + 1 



A /2m 



m 



l/2m 



1 \ i/2m 

1 + -J +o{i/j). 



and, by a straightforward computation: 

F{\Xj - EXj\ > n^) < Ce~"' 
for 1 < J < n. Consider the matrix (j3ip . with entries 

2l/2m 



aj = N{0,1). 

Let Jn = {mam/nY/^'^'^^ An, with An = EAn, and 

m 

J-i 27 

i=i ■' 



dr. 
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The matrix J„ approximates a Jacobi matrix associated with the orthogonal polynomials 
for the measure (136p . Indeed, letting 6„ be the coefficients in the three-term recurrence 
relation: 

xpn{x) = bn+lPn+l{x) + anPn{x) + hn-lPn-l{x), 

we have, by [5j Eq. (2.33)]: 

l/2m 



I /n+ 1 



2 V Om'TT. 



+ 0(1/712-1/2^). 



Note that = 0. The eigenvalues of the matrix 

^n-l hn-2 



%-2 



6l 



V 



Their 



6o 
6o / 

are the (rescaled) zeros of the nth orthogonal polynomial with respect to ([36 
limiting density can be explicitly computed; see O Eq. (2.4)]. It has the form 

/^(dx) = -x2/i„(x)l[_i^i](x)dx, 

where hm{x) is a polynomial of degree 2m — 2. A calculation as in [20] using the Riemann- 
Hilbert asymptotics in [5] shows that the normalized trace of {Jn — z)~^ approximates the 
Stieltjes transform of the limiting density fi{dx) with precision 0(l/n) for > n"^"'"^ 
and z away from ±1. By the resolvent expansion argument of Section [2] and Corollary [3l 
convergence of the empirical eigenvalue distribution for {aran)~^/'^'^An holds for intervals 
of size |/| > strictly inside the bulk of the limiting density. 
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